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SUMMARY 


Uniformly valid equations are developed for calculating spatial derivatives of 
flow quantities behind a curved shock of any strength. For steady two-dimensional 
inviscid flow, explicit formulas in terms of shock curvature are derived. Those 
factors which yield the equations indeterminate as the shock strength approaches 0 
cancel analytically, resulting in equations which are valid for shocks of all 
strengths. 

An example of an application of the method is outlined for the problem of shock 
coalescence in which asymmetric effects are included through derivatives in the cir- 
cumferential direction. Ihe solution of the coalescence problem requires values for 
spatial derivatives of the flow variables behind a resulting shock which is often 
very weak. 


INTRODUCTION 

In many flow-field computational problems, it often becomes necessary to calcu- 
late spatial derivatives of the flow quantities. Numerically determining these 
derivatives by finite difference or other approximations often leads to great inac- 
curacies in the solution. Analytical methods of calculating spatial derivatives of 
flow variables behind curved shocks are developed in references 1 to 3. These 
methods involve taking the tangential derivative of the standard shock equations 
along the shock, thus resulting in expressions relating the derivatives behind the 
shock and the curvature of the shock. This approach works well and has many applica- 
tions except when the shock strength approaches 0 and the equations become indetermi- 
nate. 


Included in this paper is a method of determining, for these derivatives, 
expressions which are shown to remain finite as the shock strength approaches 0. An 
equation for the jump in the derivative of a flow variable across the shock is 
derived by using the flow equations in section A of part I. It is the ratios in this 
derived equation which cause difficulties as the shock approaches a Mach wave. In 
section B alternate expressions, which remain finite, are derived for these ratios by 
using the standard shock equations. 

An application for the method may be found in reference 4 where the shock- 
coalescence problem with asymmetric effects is solved in a sonic boom propagation 
routine. The asymmetric effects are felt through cross derivatives of the flow vari- 
ables. It is in solving part of the shock-coalescence problem that it becomes neces- 
sary to know the derivatives of the flow variables in both the axial and the radial 
directions behind the resulting shocks. In this application, the calculations are 
performed away from the disturbing body, and the shocks involved often become so weak 
as to cause numerical difficulties. In the method developed in reference 4, the 
approximation is made that one of the resulting shocks is so weak that the flow vari- 
ables can be assumed to be continuous across the shock, thereby limiting the appli- 
cability of that method to cases in which this is true. 



The method developed in part I removes the need for a weak -shock approximation 
in the coalescence solution and thus makes the procedure applicable to shocks of any 
strength. Part II of this paper outlines the procedure for incorporating the equa- 
tions derived in this paper into the coalescence system. 

This work was done to satisfy in part the requirements for the Degree of Doctor 
of Science at George Washington University^ February 1983. 

SYMBOLS 

A,B dummy variables 

right side of equation ( 11 ) 
a speed of sound 

= a - (q • v)-^; zero for characteristic 
^ defined beneath equation (7) 

shock surfaces 

f^,f*^ curvature of shock surfaces 
Q defined beneath equation ( 8 ) 

g acceleration due to gravity 

H enthalpy 

h slipstream surface 

i unit vector in axial direction 

A 

j unit vector in radial direction 

^ defined in equation (9) 

k shock curvature 

m mass (see eq. (27)) 

n direction normal to streamline 

p pressure 

q velocity 

R gas constant 

S entropy 

f 82 ^ . . . ,S 2 q ratios defined by using shock equations (see part 1(B)) 
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T 


temperature 


U 




(j^O 


u, V 

V. 


'1 


second circumferential derivative of axial velocity in plane of 
symmetry 

velocity components in axial and radial directions, respectively 
velocity 


V 


second circumferential derivative of radial velocity in plane of 
symmetry 


4 ^ 4 ^ I 4^0 

X( (j;) ,R( ();) intersection of shock surfaces (fig. 2) 

x,r, ([; axial, radial, and circumferential coordinates, respectively 
2 vertical coordinate in Cartesian coordinate system 

a shock angle in general problem 

shock angles of f^ and f^, respectively, in coalescence problem 
ratio of specific heats 
flow angle 
arc length 

direction normal to shock 




V unit vector normal to shock 

p density 

a direction along streamline 

T direction tangent to shock 

A 

T unit vector tangent to shock 

Subscripts: 

X, r first derivative with respect to x or r, respectively 

V derivative normal to shock 

T derivative tangent to shock 

1,2 region ahead or behind shock, respectively 

4,5 regions 4 and 5 in shock -coalescence problem 



Mathematical notation: 


* first derivative with respect to x 

'* second derivative with respect to x 

vector quantity, when over symbol 
[| I ] jump in variable across shock, for example, [|^| ] “ ^2 ** 

*> approaches 


PART I 

DERIVATION OF EQUATIONS FOR SPATIAL DERIVATIVES BEHIND 
SHOCK IN TERMS OF CHANGES IN SHOCK ANGLE 

A - Development of Ratios Needed Across Shock and Explicit 
Formulas for Derivatives 

There are essentially five steps involved in determining the formulas for the 
needed derivatives. For ease in following the procedure, the steps are outlined 
next. After the outline, the details of the derivation are given in each step. 

(1) Resolve the flow equations into derivatives parallel to, 5/5 t, and perpen- 
dicular to, 5/5v, the shock. 

(2) Reduce these equations to one equation in terms of 5/5v of one variable, 
say u^; that is, 


C (5u^/5v) = R.H.S. 
N^ 2 


where R.H.S. represents the right-hand side of the equation. 

(3) Repeat steps (1) and (2) for the field ahead of the shock. Form a differ- 
ence between the two equations: 


c^J|au/av| ] + (&a^/5v)[|c^| ] = [Ir.h.s. I ] 

where [| | ] represents the jump across the shock. 

(4) Solve for [] 5u/5v| ]; that is, 

[l5u/avl] = (-au^/av)[|c^iyc^^ + [\R.K.s,\]jc^^ 
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The shock equations are used (see section B of part I) to show that the numera- 
tors of the ratios and [ |r.H.S. | are divisible by Cj, . Conse- 

quently, explicit formulas can be derived for those quotients which will be finite as 
the shock strength and therefore approach 0. 

(5) Derive explicit formulas for the derivatives behind the shock and the curva- 
ture of the streamline in terms of the slope and the rate of slope change of the 
shock and of the local flow field and the derivative ahead of the shock. 

Step 1 ♦ - As shown in figure 1 , a is defined to be the shock angle with respect 
to the x-axis and 0^ is the flow direction^ where the subscript j = 1 indicates 
flow ahead of the shock and j = 2 indicates flow behind the shock. (The num- 
bers (1) to (5) in the figures refer to regions 1 to 5.) 



a>0 


a<0 


Figure 1.- Coordinate system. 


With these definitions, the unit vectors tangent and normal to the shock become, 
respectively. 


A ^ /V 

X = i cos a + j sin a 


( 1 ) 


A A A 

V = (i sin oc - j cos a) sgn a 


( 2 ) 


(where sgn represents the sign (signum) of the argument) and the first partial 
derivatives with respect to x and r in terms of x and v are, respectively. 


5x 


5x 


cos 


S 

a + sin a sgn a 


(3) 
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5 5.5 I 

sin a - cos a sgn a * 

5r 5x 5v 

The flow equations given in cylindrical coordinates x,r (see ref. 5) in the 
vertical plane of symmetry and rederived in reference 4 are 


2 

(a^ - u^)u - uv(v + u ) + (a^ - v^ )v - - - gv 

X X r r r 


V - u = 
X r 


1 5H ^ T 5S 
q 5n q 5n 


where q(5/5n) = u(5/5r) - v(5/5x), and 5/5n is the derivative normal to the 
streamline. Applying equations (3) and (4) to equation (5) and rearranging yields 

[(a^ - u^) sin a + uv cos a]u sgn a - [uv sin a + (a^ - v^) cos a]v sgn a = 


where 


^ + gv + r (a^ - u^) cos a - uv sin alu + F-uv cos a + (a^ - v^) sin alv 

r L j L j ^ 

Likewise after applying equations (3) and (4) to equation (6) and rearranging, one 
obtains 


(u cos a + V sin a) sgn a = -Q 
V V ^ 


where 


9 4 ( 11 -^ 11 )-" 


Sin a + V cos a 
T X 


Step 2.- Elimination of v from equations (7) and (8) yields the following 
single equation for u^: 


[a^ - (u sin a - v cos sgn a 
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or 


[a2 - (5 • = -9< sgn a 

where 


^ ^ sin a + Q [uv sin a + (a^ - v^) cos a] 

As the shock approaches the characteristic, a^ approaches (q • v)^ and hence 
the right-hand side ^ which goes to 0, becomes the characteristic equation. This 
explains why the equations for the derivatives become indeterminate as the shock 
strength approaches 0. 

Step 3 .- Equation (9) applies to the regions ahead of and behind the shock. 
Applying it to each and forming their difference gives 


A ^ 2I ^^^2 
5v 


If [| |] = Jump, for example, [|*3<|] “*^ 1 ' equation (10) may be 

rewritten as ^ ' 



v)2| ] = -[l£Kl ] sgn a 


where use has been made of the fact that 

[l“l] - ■ ’'2°2 - *2®1 ^ ■ ^[IbI] + Bjlal] 

Step 4» - By solving for [l5u/5v|], one obtains 
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It is now necessary to show that the right-hand side of equation (11), or 

remains finite as a^ - (q • v)^ ^ 0. If is defined as a^ - (q • v)^, 

2 2 ^2 2 2 

the first term of the right-hand side of equation (11) becomes 


1 - 


- (q^ • 0)2 


N. 


To compute [|‘tK|]y^Cjj , expand [1*0(1] to get 


[ 1 0< I ] = [ I I ] sin a + [ I 9 I “] 


2 2 ' 


+ 9^([|uv|] sin a + [ja^ - | ] cos a) 


Expanding further gives 


n ^l] 1 [|a^v|] _ g[|v|] l[\a^ - u^|] 

C. 


[h'l] . 


-XT r C C 

2 2 2 


COS a sin alu^ 


r 2 2 

[ (a^ - ) cos a - sin ccj-^ 


_/[M] _ „ _ [l 


2 2 n 

a - V I J 


cos a - 
^2 **2 


sin a\v. 


- [v^u^ cos a - (a^ - v^) sin a] 




and 


[191] _ [|(1/q)/(5H/5n)|] [|T(5S/5n)|] 


[Kt;|] 

- sin a — + cos a — 


N, 


N. 


N, 


N, 


then 


( 12 ) 


(13) 


(14) 


(15) 
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Therefore, to evaluate , it is necessary to evaluate and show that 

numerous ratios remain finite as 0. These evaluations are provided later by 

the shock equations in part 1(B) of this section as indicated by the notation 
beside the ratio. These ratios are 


[|v|] 


C " 
^2 


[|uv|] [l^t] [Kl] 

±J LL = iJ Li V + — Li u = S V + S u 

C C 2 C 1 6 2 7 1 

N N N 

^2 2 2 


[|a' - v^l] _ . 

C 14 

^2 


[la" -u^ll _ c 

C ®15 

^2 


[|a"v|] _ s 

^^2 


[|(1/q)(5H/5n)|; 


[ |(T/q)(dS/an)l] 


n \ i ] 

— = s s 
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To evaluate expression (12), the first term on the right-hand side of equation (11), 

2 -,2 

a. - (q^ • v) 

is shown to remain finite by S^. 

^2 

With the definitions of the aforementioned ratios and the determination that 
they remain finite, then the right-hand side of equation (11), is now defined 

and is finite for all shock strengths • 

Step 5 .- The explicit formulas for the derivatives in region 2 behind the shock 
are now developed. All flow quantities and derivatives in region 1 are known and all 
flow quantities (not derivatives) in region 2 are known. 

From equation (43), in the equations developed after this section. 


^ • T = S 2 k[a^ " (q 2 * sgn a (16) 

and from equation (46), 

Jlf|_ • '■ ■ ■ ‘^2 • 

By definition. 



Also by definition. 
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By equating the i and 3 components in the aforementioned definitions, the follow- 
ing relationships are obtained. From the 1 component, 


hi. 




5q 


• T COS a + 


8 q 

5t 


• V Sin a 


bn 


2 


bx 



bx 


X cos a + 


5t 


V sin a 


By using equations (16) and (17), this becomes 


bn^ 3u^ 

— = -r — + S^kC„ cos a sgn a + C„ sin a 
bx bx ^ ^2 ^2 


(18) 


Likewise from the j component. 


a[|v|] 

1 5q 

bx 

J ax 


X sin a - 


5q 
bx i 


• V cos a 


and thus 


9v2 5v^ 

■r — = T — + S^kC sin a sgn a - C cos a 
bx Qt 2 ^ 11 

2 2 


(19) 


From equation (11), one obtains 


au2 au^ 


+ lA 


bv bv 
and from equation (8), 


( 20 ) 


^ 

bv 


sin 


a 02 


(21 ) 


By using the derivatives in terms of x and v, the following relationships may be 
defined directly: 


^^2 ^"2 

5x bx 


cos a + 


5u. 

bv 


sin a sgn a 


( 22 ) 
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sin a 


cos a sgn a 


(23) 


5r 9 t 

!l2 ^ ^ 

9x bi 
5r 5T 


cos 


sin 


9u, 

z 

5v 


^'’2 . 

a + sin a sgn a 


a + cos a sgn a 


where these four derivatives are functions of the shock curvature k. 


(24) 


(25) 


B - Definitions of Ratios to S 2 Q by Using Shock Equations 

The jump in the normal derivative across the shock is given by equation (11). 

The evaluation of this equation is defined in terms of ratios to ^ 20 * In this 

section the standard shock equations will be used to define the ratios and to show 
that they remain finite as the shock approaches a characteristic. 

By rewriting the shock conditions as derived in reference 2 in terms of the jump 
notation [ | | ] , the r momentum equation becomes 


[|q*T|]=0 (26) 

the continuity equation becomes 

[|pq • v|] = 0 . 0 = m = p^q 2 * V (27) 

the X momentum equation becomes 

[|p(q • V)^|] + [1p|] = 0 m[|q • v|] + [jpl] = 0 (28) 


and the energy equation becomes 


[|(q • v)^|] + - "1 - ^ r[|t|] = 0 


By using the equation of state 


(29) 


p = pRT 


(30) 
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and equations (27) and (28) to eliminate p and p, the following equation is 
obtained: 


[|q • v|] + R 


= 0 


q • V 


Expanding [ |T/(q • v)|] and rearranging gives 



T 



— ^ 
q • V 



~ 


^^2 * ^2 ' ^ ^2 ' ^ • V '^2 * ^ ^^2 * 




Therefore, equation (31) may be written as 

[|q * v|][(q^ • ^^^^2 • v) - RT^] + R[|T|](q^ . 0) = 0 

Equation (29) may be rewritten as 


(q^ • V + q^ • v)[ |q • v| ] + R[ |t| ] = 0 

Eliminating [|t|] from equations (32) and (33) gives 


[|q • v|] 


and since [ |q • v|] ^ 0, 


T (q^ • V)(q 2 • v) - ^ RTi - (q^ • v)^ - (q^ • v)(q 2 • v) 


= C 


^ mm yV 


2YRT-I = (y + 1)(q^ • v)(q 2 • v) - (q^ • v) (y - 1 ) 


or 


(y + 1)(q^ • v)[jq • vj] = 2[yRT^ 


- (q^ • v)^] 


which becomes 


(Y + inq^ • v)[|q • v|] = 2a J 


1 - 


(q^ * v)‘ 


(31 ) 


0|] 


(32) 


(33) 


(34) 


(35) 
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or 


[ |q • v|] 2 ^ 

- (q^ • v) (y + 1 ) (q^ • v) 

Equation (32) may be rewritten with 

Hi * v|]t2 ^ [|t|] 

(q^ * ^^^^2 * q^ • V 

to obtain 


2 yRT 2 = 2(q^ • v)^ - (y + l)[lq * v|](q 2 * v) 
or 


[ |q • v|] ^ -2 ^ g 

^2 ■ ^^2 * ^ 

From equations (36) and (37), one obtains 




From equations (33) and (37), one may obtain 



and from equations (28) and (37), 

[ |p|] ^ 2m ^ ^P2 

a^ - (q2 • v)^ (y + ^)(q2 * 


(36) 


(37) 


(38) 


(39) 


(40) 
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Note that 


[|5 • t|] = [|u|] COS a + [|v|] sin a = 0 

[15 • v|] = ([|u|] sin a - [|v|] cos a) sgn a 

[ I'^l ] = [ k * ■’^1 ] COS a sgn a + [k • v|] sin a 

[k|] = [|q * T^l] sin a - [k • v|] COS a sgn a 


where the first expressions defined in [|u|] ®"d [|v|] go to 0. Therefore, from 

equation ( 37 ) , 


[k|] 


2 2 


- = sin a = Sg 


(41 ) 


[|v|] 


• ^^^2 * 


= -S„ COS a sgn a = S_ 
^.2 2 7 


(42) 


The ratios to Sj are all used for the v derivatives behind the shock. 

From the aforementioned expressions, it is seen that these ratios all remain finite 

as the shock strength approaches 0 ^i.e., a^ - (q 2 * -»• o) . 

The tangential derivatives of [| j] across the shock will now be developed. 
If T is defined as the arc length along the shock, then, from equation (26), 


• T + kv[|q| 


since 


fetii 


n - 0 - [Ilf ’ll • ^ B ■ 


Qq 

5t 


y\ 

dt 


A 

+kv 


A, 

&v 

5t 


A 

±kx 


where k is the curvature of the shock. Therefore, 



• X = ±[|q • v|]k 


15 



and 



&q 

• T 


dx 

• T# 


a - (q • v) 
2 2 


2 = ±52^ sgn a = Sg 


(43) 


From equation (37), 


fe[l5-l] = 


9q 


5 t 


V ± [ |q| ] • kT 


and since [|q|] • T = 0/ then 


=l] - 


5q 

5t 


Taking the 5/5 t of equation (36) yields 


Sq 


V = - 


1 


a(q^ • v) 


Y + 1 /“ ^v2 

* (q^ • v) 


9t 


^ 


v)‘ 


2 1 5 r 2 ^,2 

— J 


• 2 — ^ 2 
and dividing by a^ - (q • v) gives 



&q 

A 


_ 

. Sx . 

J * ^ 2 

1 


d(q^ • v) 


2 ^,2 
- (q^ * v) 


Y + “I /“ ^\2 

' (q^ • v) 




Y + 1 5 . 0 2 ''.2 

^1 V - (q^ • v) 


= S, 


(44) 
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To show that 


a 

5t 


a 


2 

1 


a 


2 

1 







remains finite as a^ - (q^ • v)^ -> 0/ the following substitutions will be made. 

Denote \ as the arc length along the characteristic, in the same sense of z 
along the shock, and as the angle of the characteristic with the x-axis . Then, 




= 0 


a a a 

^ = cos ai ^ + sin «i ^ 


and 


ar 


cos a — + sin a — 

ax ar 


a ^ 

ax " a\i 


(cos a - cos )^ + (sin a - sin 


2 — 2 

Recall that as a ^ • v) 0. The equation 


cos a - cos a, = 


1 2 2 

. ; (cos a - cos a.) 

1 cos a + cos 1 


can be rewritten as 


cos a - cos = 

1 cos a + cos a. 


2 2 
-(sin - sin a) 
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Multiplying and dividing by gives 


cos a - cos = — 

(cos a + cos ) 


2 

|a, - (q, 


since sin and q sin a = q • v. Similarly, 

1 1 1 ^1 


2/ . 
(si 


: — ' — ^ — r[‘^i 

in a + sin 


2 

v) - a^ 


Therefore, 


2 ''.2 '>.2 2 
a 5 °i - ’■^1 • 5 

6T " 2. . dx ^ 2, . , . 5r 

q^ (cos a + cos a-| ) q^ (sin a + sin a-| ) 

By substituting this value into the second expression of equation (44), the following 
results are obtained: 


-k - <5, • ‘'”'1 - '5, • v>'l - — k - (5, 

dTl^ 1 J axL 1 1 J a\A 1 i 


- (q^ • v)2 


2 -.2 

- (q, • v) 


a fz '‘. 2I a r 2 


q^ (cos a + cos a^) 


q^ (sin a + sin ) 


Then, S^q remains finite as a • The term 5(q • v)/9x in equation (44) may 

be expressed as 


5(q * v) _ 5(q • v) 
5x " dx 


5(q • v) 


cos a + ^ sin a 
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where 


5 (q • v) _ 5(q • v) ^ f 5q 
5x 


5v 


dT ^ = FT * ^ ^ “ 


3q - 

• V + sgn a Oct • q) 


sec a 


and 


5(q ^ v) 
5r 


/ Bq 5v 

= a? * 


CSC a = 


Bq ^ A — 

• V + sgn a (kx • q) 


CSC a 


Therefore, 


B(q • v) J Bq 


+ sgn a (kx • q) 


Equation (44) becomes 


[ ISq/Stl ] » V _ 4 ^ 


,2 -,2 

- (q^ • V) 


Y + 1 /- ^.2 

(q^ • v) 


'5q^ 

idx 


• vj + (kq^ • x) sgn a 


+ 


^ 

(y + 1 )(q^ 



'11 


(46) 


where S^-j remains finite. 

By definition, 

Bq ^ 

• X = cos a + sin a 

Bq ^ . 

• V = (Uq; sin a - cos a) sgn a 
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Thus, 


I 6q 

Ut- = COS 


- (i • =) 


sin a sgn a 


and 


/ 5 q ^ \ . / 6q A 

~ • '^j ^ ” 1“^ • cos a sgn a 


[ 1^1:1] [iQq/a-rl] • T [|aq/aT|] • v 


COS a + 


sin a sgn a 




N. 


= S cos a + S 
8 11 


sin a sgn a = S 


12 


Similarly, 


[|vj] 


—Q = Sq sin a - cos a sgn oc = S-i^ 


[|a I] - af yr[|t|] 


^N 2 ^N 2 ''N 2 


C„ “ 


[|v2|] (Vi + V2)[|v|] 


'N- 


= (v^ + V2>S7 


[|a2 - v2|] [|a2|] - [jv^l] 






= YS4 - (v^ - V 2 )Sj = S 


14 


[ - u2|] 


-C- YS 4 - (Ui + U2>Sg = 3^5 


and 

[|a^v|] [|a^|]v2 2 

cT " c!| ‘'^®4'^2 ■'■ " ®16 

N_ 


(47) 


(48) 


(49) 


(50) 


(51 ) 
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By definition 


[|h|] = 0 




— - &H/5t 

9n sin(a - 0) 


in regions 1 and 2, and 


5n 


sin(a - 0) 


= 0 


sin(a - 0) = (q • v)/q 


Therefore, 


1 9h,- 

i aTT*'’ • '■> 


= 0 


Expanding the previous equation gives 


1 M 

q 5n 


• ^>2 q 


HS),tli-0|]-o 


When dividing by Cjj / one obtains 



1_ ^ 

1 


1 

. 

q 5n 

; 

dn [|q • v|] 

dn 


(?2 • 

v)2 ^2 


?2 • 0 


"2 = "17 


(52) 


Also by definition in regions 1 and 2 , 


5 S 

9a 


0 
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and 


5S 5S/5t 


q 5S 


5n sin (a - 0) q • v 


Thus, 


1 as,- as 

i 5^' sT 


T ^ 


T(as/ax) 

q an 


q • V 


Expanding t by using the expression for entropy in reference 6 gives 
or 


5S 


^ R 

aT 





5t 

■ 

Y - 1 

St 

- 

p ax 

- 


T 

5S 

_ Y 


1 

aT 

1 

5p 

A 

• V 

St 

Y - 

1 ^ 

q • 

0 at 

Pi (5i * V) . 

St 


After dividing by , 


1 5T 

5t 


q • V 


f C and 
^2 


ap 


St 


'c must be evaluated, 
^2 


To examine the factor [ 1 I / recall from equation (28) that 

[1p 1] + 41^ • v|] = 0 


Taking the tangential derivative and dividing by gives 

^2 

[|p|]^ + \[|i * v|] + m[|i • v]]^ = 0 


[|p^l] [|q • v|] m[|q ♦ v|]. 


= -m, 




C. 


N. 


N. 


= -m 


2 ‘*2 

[15 • v\] [lil]^ • V 


c 


- m 


N. 


N. 
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Thus 


[ipj] 


N. 


= -m 

T 2 11 18 


(53) 


The factor 


1 5T 


|q • V 6x 

r[|t|] from the energy equation (eq. (29)): 


shall be examined by first obtaining an expression for 


^[|t|]^ = - • v)^l] 


Y-1 


[ I (q • V) (q • v) _^1 ] 


= - Y - I 


(q 2 * v)^[|il]^ • vj + (q^ • v)^[|q • v|] 


After dividing by , this becomes 


r[|t|] 


^ = ♦ (5, • v),S2] 


By applying the expansion formula for []ABj] (shown below eq. (10)) to 
r[ |t /(q • v) |], one obtains 


R 


q • V 


r[|t^|] 
^2 * ^ 


+ RT 


q • V 


-[|5-'v|] 


+ RT 


^2 • ^ 


1 


(q.| • 
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Therefore, after substituting for r[|t ] ]/(q • v), 

T 2 


r[ |t /(q *0)1] 


Y - I 


11 




RT 


1 


• v)(q. 


V) ^2 


= . Y - 1 


(q • v)(q * v) - RT, 
1 1x1 


Sl1 + ^2 


(q^ • vMq^ • V) 


to the 

evaluation of 

T ds 

. q 

r 

3S /,- 

c 


CN 


= - 

■;=ta 


[\V4K 


gives 


•0 


9p 


5t 


Y - 1 


RT^/(q 


N, 


= S 


20 


= S 


19 


;>|] 


( 54 ) 


(55) 


The quantities S-| to ^20 finite as approaches 0 and thus provide 

uniformly valid expressions for the ratios needed to evaluate equation (11) and 
therefore 5u /5v. The solution process for the derivatives behind the shock con- 
tinues with step 5 where the spatial derivatives are determined as functions of the 
curvature of the shock. 


An application for this system of equations will be shown in part II, in which 
the equations for shock coalescence with asymmetric effects are developed. 


PART II 

APPLICATION OF METHOD TO SHOCK COALESCENCE WITH ASYMMETRIC EFFECTS 

The coalescence system is briefly described here. (See ref. 2 for a detailed 
derivation of the governing equations of the system.) Two shock surfaces f (r,4;) 
and f^(r,(\)) coalesce to form a resultant shock surface f^(r,(];), a contact^ surf ace 
h(x,cp), and a weak (isentropic) shock or expansion or the opposite family f (r,4>). 
(See fig. 2.) The problem is first solved axisymmetrically, and at some distance 
away from the body it may be treated as a two-dimensional problem, as represented in 
figure 3. By assuming that all conditions in regions 1, 2, and 3 are known and that 
the point of coalescence is known, the unknowns for the system are u^, v^, p^, 

^ 5 f v^, p^, T^, and where p is the shock angle with respect to 

the upstream flow. The equations for the system are four shock equations across f^, 
four shock equations across f^ (assuming that f^ is a shock), and matching pres- 
sure and flow direction at the slipstream surface h. 
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Resultant shock 



Figure 2,- Three-dimensional shock coalescence* 


f^(r ) - Resultant shock 



Figure 3*- Two-dimensional representation of shock coalescence. 


The system is closed with 10 unknowns and 10 equations. To derive the asym- 
metric system" of governing equations, the system is treated as five intersecting 
surfaces along X{(\>) ^ (See fig. 2.) The governing equations are the five 

intersection equations, five shock equations at f^ (continuity, energy, and three 
momentum), five shock equations at f^, and three slipstream equations at h. The 
second derivative with respect to <\> (circumferential direction) is taken for each 
of the governing equations along the corresponding surface. The equations are then 
reduced to the (j; = 0 plane resulting in 18 unknowns of the form ^44|(p=0' 

and so forth. The asymmetric system would be closed with 18 unknowns and 18 equa- 
tions, except that it becomes necessary in deriving the asymmetric equations to know 
the spatial derivatives of u, v, p, and T in both regions 4 and 5. 
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The method derived in part I of this paper is applied and equations (22) to (25) 

provide a method for obtaining u. , v. , f r ^^id v^ in 

^r ^r ^x ^x ^r ^r 

terms of f^ and f^, respectively, which are unknown. 

The Euler equations in the plane of symmetry, 4; = 0, are valid in regions 4 
and 5 and are now applied to find relations for p^^ and in those regions. 

Using the Euler equations for p^ and p^ gives 


P 4 = 

X \ X r 


P5 = ~/-P' 

X 


5U5U5 - P3V5U5 \ 

X r i 


-p.u.v. - p.v.v. - 

4 4 4 4 4 4 


r T. \ X 

4 




P 5 - ^('-P5»5''5 - Ps''5''5 - 9 P 5 

r T^ y X r 


(56) 


(57) 


(58) 


(59) 


The three equations valid at the slipstream surface h (the tangential derivative 
along h of slipstream conditions) are 


p + p h = Pc + Pc ^ 
^4 4 X 5 ^5 X 

X r X r 


(60) 


u^ h + u^h 
4 X 4 XX 
X 


v^ + u^ h - V, h 
4 4 X 4 X 

X r r 


= 0 


(61 ) 


Up h + u-h 
5 X 5 XX 
X 


v_ + Up h - V- h 
5 5 X 5 X 

X r r 


= 0 


(62) 


There are 15 equations which may be solved for the 15 unknowns , u. , , 

X r X 

'^^x' '^^x' ^'^x' ^^x' ^3' ^4' 

To determine the derivatives for temperature T, the enthalpy relationships are 
used, where enthalpy H is defined as 


H 


X 


Y " 1 


RT + 


gz 
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By definition, 5H/9a = 0 where 




and from the shock equations derived in the previous section, 

[liiO 


= s 


17 


Therefore , 


/ , an ' 

= q, c„ s,, + — 


5n 


^2\ 17 5n 


From the definition of enthalpy. 


5n Y - 1 2 5n 


dT„ 




Y - 1 5n 


+ u 




2 3n ^2 5n 


Thus, 


5n Y “ *^2 5n ” '^2 9n 


V- ira 
Y \dn 


"2 ^“2 " 2 ''2 °“2 

5x 


5r dx 


Along the streamline, 


, / au 3v 

2 Y - 1 2 2 

5a Y V ^2 ao ’^2 ao \ 


= Y - 1 ( ^ ^ ^ 

Y \ q2 q2 ^2 ^2 


(63) 


(64) 
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From the normal and streamline derivatives of T, the following relationships are 
determined: 


(65) 


( 66 ) 


Equations (65) and ( 66 ) applied in regions 4 and 5 provide the four additional equa- 
tions needed to determine T* , » Tc » and Tr • 

^r X r 

For the curvature of the streamline ^ ” q^(x) behind the shock, note 

that ^ 



^2 

*>■'2 ''2 

dT2 

Sx 


&a ■ q^ 

Sn 

9T^ 

v^ 

ST- u^ 

ST- 

2 

2 

2.2 

2 

Sr 

^ ^2 

Sr q^ 

Sn 




^ {1 + 


where ^ ^2 ^ ^2 

streamline r = 92 (x) yields 


The derivative of this equation along the 


7 + V q* = U q" + 

2 2 ^2 2^2 
X r 


u 


+ u q * W * 
2 2 ^2T2 

X r 


Therefore, 


k 


2 



+ 




2 / 2 ) «2 


u|v2^ <■ 'l;V2(v2^ - - <12/2 


(»2 ^ v 2 ) V2 


(67) 


All quantities on the R.H.S. of equation (67) are known behind the shock from equa- 
tions (23) to (25) and (56) to (62). Thus, the curvature of a given streamline 
behind the shock may be found by using equation (67), 
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CONCLUDING REMARKS 


A set of uniformly valid explicit formulas for spatial derivatives of flow quan- 
tities behind a shock wave have been developed in terms of the shock curvature* The 
equations for jump conditions across a shock have been rearranged so that those quan- 
tities which become indeterminate as the shock strength approaches 0 have been elimi- 
nated analytically. The resulting formulas are therefore valid for shocks of any 
strength. 

In the solution of shock coalescence with asymmetric effects as applied to sonic 
boom propagation/ the need for spatial derivatives behind the shocks arises. Shock 
strengths in sonic boom problems may vary from moderate to very weak. This problem 
then has need of equations which are valid for finite shocks and also for shocks 
approaching zero strength. An outline is provided to show the application of the 
uniformly valid formulas developed in this paper to the shock-coalescence problem. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 21, 1984 
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